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ABSTRACT 

Variational  bounds,  both  upper  and  lower,  are  found  for  the  equivalent  spring 
constant  of  a  double-strap  joint  which  represents  a  sub-element  of  bonded  repairs  to 
cracked  structure.  Conservative  estimates  of  the  equivalent  spring  constant,  needed  for 
accurate  design,  are  obtained  from  variational  analyses  of  the  joint.  Estimates  from 
various  analytical  models  of  varying  level  of  approximation  were  obtained.  Simpler 
expressions  for  the  spring  constant  resulted  from  relaxing  certain  assumptions, 
however,  the  theoretical  guarantee  of  a  true  upper  or  lower  bound  was  lost.  Spring 
constant  estimates  were  compared  with  finite-element  model  results  and  so  the  fidelity 
of  the  variational  bounds,  specially  for  the  simplified  analyses,  could  be  established. 
An  improved  formula  is  proposed  for  use  in  design  procedures  in  RAAF  C5033. 
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Variational  Bounds  for  the  Equivalent  Spring 
Constants  for  Bonded  Repairs 


Executive  Summary 


Bonded  composite  repairs  are  an  effective  and  cost-efficient  means  of  restoring  the  life 
of  cracked  metallic  aircraft  components.  The  Aeronautical  and  Maritime  Research 
Laboratory  (AMRL)  has  used  this  technique  for  over  25  years  to  repair  damaged  Royal 
Australian  Air  Force  (RAAF)  aircraft.  A  critical  step  in  the  repair  design  process  is  to 
obtain  an  accurate  estimate  of  the  stress-intensity  factor  associated  with  the  crack  in 
the  metallic  component  after  application  of  the  repair.  In  the  initial  design  stage, 
analytical  formulae  are  used  to  obtain  these  estimates. 

This  report  evaluates  various  analytical  models  (some  developed  by  the  authors)  that 
can  be  used  to  obtain  post-repair  stress-intensity-factor  estimates  through  the  concept 
of  equivalent  spring  constants.  The  bonded  composite  repair  can  be  thought  of  as  a 
series  of  springs  resisting  the  opening  of  the  underlying  crack  when  a  load  is  applied. 
The  springs  are  characterised  by  a  spring  constant  and  the  validity  of  the  various 
models  is  explored  through  consideration  of  this  constant.  Of  particular  concern  is  that 
the  models  provide  truly  conservative  estimates  (an  upper  bound  for  the  spring 
constant)  that  can  reliably  be  used  in  the  design  process.  For  this  reason,  a 
complementary  energy  functional  for  the  bonded  repair  was  constructed  with  the  use 
of  either  stress  or  displacement  based  potentials  and  a  variational  analysis  undertaken 
to  provide  an  upper  or  lower  bound  respectively. 

The  most  suitable  model  for  design  purposes  -  one  which  provides  conservative 
design  estimates  and  is  not  too  mathematically  complex  -  was  identified.  This  model 
will  be  proposed  to  the  RAAF  for  inclusion  in  RAAF  Engineering  Standard  C5033  - 
"Composite  Materials  and  Adhesive  Bonded  Repairs". 
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1.  Introduction 

Bonded  composite  repairs  are  an  effective  and  cost-efficient  means  of  restoring  the  life 
of  cracked  metallic  aircraft  components.  The  Aeronautical  and  Maritime  Research 
Laboratory  (AMRL)  has  used  this  technique  [1]  for  over  25  years  to  repair  damaged 
Royal  Australian  Air  Force  (RAAF)  aircraft.  A  critical  step  in  the  repair  design  process 
is  to  obtain  an  accurate  estimate  of  the  stress  intensity  factor  associated  with  the  crack 
in  the  metallic  component  after  application  of  the  repair.  In  the  initial  design  stage, 
analytical  formulae  are  used  to  obtain  these  estimates. 

The  purpose  of  this  paper  is  to  assess  the  accuracy  of  the  current  technique  and  to 
explore  other  analytical  techniques  that  might  provide  improved  estimates.  Truly 
conservative  estimates  are  desirable  for  the  initial  design  stage  and  the  conservative 
nature  of  the  estimates  is  assessed  through  the  concept  of  variational  bounds. 
Variational  bounds  for  the  equivalent  spring  constant  of  a  bonded  repair  arise  from 
minimising  the  complementary  energy  functional  when  either  stress  (lower  bound)  or 
displacement  potentials  (upper  bound)  are  employed. 

The  accuracy  of  the  various  models  is  established  by  comparison  with  finite-element 
analysis  results  that  are  taken  as  "exact".  A  second  criterion  by  which  the  models  are 
judged  is  the  simplicity  of  the  expressions  obtained  for  the  design  parameters  of 
interest.  Solutions  based  on  prescribed  stress  potentials  that  satisfy  the  equilibrium 
equations  of  continuum  mechanics,  for  instance,  may  yield  accurate  values  for  the 
equivalent  spring  constant  of  the  repair  scheme  but  may  be  very  complex  to  evaluate. 
By  relaxing  certain  assumptions  simpler  but  still  accurate  solutions  may  be  obtained. 
However,  the  theoretical  guarantee  of  a  true  lower  bound  is  lost. 

An  important  aspect  of  establishing  the  accuracy  of  the  various  analytical  models  was 
to  consider  isotropic  and  orthotropic  adherend  joints  having  a  range  of  adherend 
thicknesses.  Some  models  are  accurate  for  only  a  small  range  of  adherend  thicknesses 
so  it  was  necessary  to  consider  a  wide  range  adherend  thicknesses. 

The  plan  of  the  report  is  as  follows:  in  Section  2  the  double-strap  joint  (Figures  1  and  2), 
which  is  representative  of  typical  bonded  repairs,  is  described,  the  concept  of  the 
equivalent  spring  constant  is  summarised  and  the  existence  of  variational  bounds 
established;  in  Section  3  the  finite-element  model  and  results  for  both  the  isotropic  and 
orthotropic  adherend  double-strap  joints  are  tabled  and  in  Section  4  the  various 
analytical  models  and  results  are  given  along  with  comparison  of  the  finite-element 
results.  Conclusions  are  given  in  Section  5. 
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2.  The  Equivalent  Spring-Constant  Concept 


2.1  Geometry  of  the  Double-Strap  Joint 


The  double-strap  joint  (DSJ)  represents  a  section  through  the  central  section  of  a 
double-sided  repair  as  shown  in  Figure  1.  It  is  also  representative  of  one-sided  repairs 
in  which  there  is  substantial  restraint  against  out-of-plane  bending. 


Figure  1.  The  DSJ  as  a  Subelement  of  a  Double-Sided  Repair. 

The  nomenclature  for  the  joint  geometry  is  shown  in  Figure  2  below: 


Figure  2.  Nomenclature  for  the  DS  J. 
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The  subscripts  i,  o  and  a  refer  to  the  inner  adherend,  the  outer  adherend  and  the 
adhesive  respectively.  The  material  properties  are  similarly  designated  with  Eir  E0  and 
Ea  representing  the  Young's  moduli  for  example. 

For  simplicity  of  mathematical  modelling  an  infinite  (L->  oo)  DSJ  was  modelled.  DSJ's 
of  sufficiently  long  overlap  length  are  well  approximated  by  such  a  model.  A  quarter 
model  of  the  infinite  DSJ,  taking  account  of  the  symmetry  of  the  joint  across  its  mid¬ 
plane  is  shown  in  Figure  3. 


Figure  3.  Quarter  model  of  the  identical  adherend  infinite  DSJ  (NB.  The  z  and  x  axes  are 
axes  of  symmetry). 

2.2  Definition  of  the  Equivalent  Spring  Constant 

Rose  [2]  introduced  the  concept  of  an  equivalent  spring  constant  to  model  partial 
reinforcement  of  cracked  substrates.  The  constraint  imposed  by  the  bonded  composite 
repair  on  the  opening  of  the  crack  is  simulated  by  a  distribution  of  springs  across  the 
crack  faces.  A  simplified  assessment  of  the  repair  efficiency  arises  from  such  an 
approach.  The  springs  across  the  crack  faces  have  an  equivalent  spring  constant,  k, 
which  is  given  by: 


cr.  -  E.ku  , 

0  i 


(1) 


where  cr^is  the  applied  stress,  E,  is  the  Young's  modulus  of  the  inner  adherend  and  u  is 
the  average  displacement  of  the  inner  adherend  in  the  x  direction. 

Equivalently,  k  can  be  defined  in  terms  of  the  complementary  energy,  Uc  [2]: 


tt  1  » 

U  =—  ( TnU 

c  2  0 


_la0 


2  E.k 

i 


(2) 
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The  complementary  energy,  when  the  strains  are  expressed  in  terms  of  stresses,  is  also 
given  by: 


(3) 


where  the  summations  are  over  all  the  stresses  present  in  the  inner  adherend,  outer 
adherend  and  adhesive  respectively. 

2.3  Variational  bounds 

As  discussed  in  Hashin  [3],  a  Uc  functional  associated  with  admissible  stresses  leads  to 
a  lower  bound  on  die  longitudinal  Young's  modulus  for  a  cracked  laminate.  Similarly 
for  the  DSJ,  a  Uc  functional  associated  with  admissible  stresses  leads  to  a  lower  bound 
on  k  while  a  Uc  functional  associated  with  displacements  leads  to  an  upper  bound  k: 

k  <k<k  (4) 

a  u  w 


3.  Finite  Element  Results 

The  finite  element  (FE)  package  PAFEC  [4]  was  used  to  model  the  DSJ.  A  plane-strain 
two-dimensional  model  (Figure  4)  employing  eight-noded  isoparametric  elements  was 
used. 


Z 


<■ 


L 


a 


u 


X 


Figure  4.  FE  quarter-model  of  the  DSJ  showing  loading  and  constraints. 
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The  following  linear-elastic  material  properties  were  used: 
adhesive 

Ga  =  700  MPa 
=  0.35 

isotropic  adherends 

£  =  71  GPa 
v  =  0.33 

orthotropic  adherends 

El  =  208.3  GPa 
Et  =  25.4  GPa 
Glt  =  7.24  GPa 
vir=  0.183 
vTL  -  0.1667 


The  orthotropic  material  properties  of  a  typical  fibre-composite  material  have  been 
simplified  to  those  of  a  transversely  isotropic  material  to  simplify  the  analysis.  The 
isotropic  plane  is  taken  to  be  the  plane  perpendicular  to  the  fibre  direction  and  only  a 
unidirectional  layup  has  been  considered. 


The  following  boundary  conditions  (see  Figure  3.1)  were  applied  to  the  model.  A 
uniform  stress,  ou,  was  applied  at  the  end  of  the  joint  and  the  longitudinal 
displacements  there  constrained  to  be  uniform  across  the  thickness: 


a  \x  =  L,0  <  z  <h.  +  h  +h  =cr 
xx '  i  a  o'  u 


(5) 


u  Mx  =  L,0  <  z  <  h.  )=u  a\x~  L,h.  <z<h.+h 
x  v  i>  x  '  i  i  a 

=  u  °\x  =  L,h.  +  h  <z<h.  +  h  +h 
x  x  i  a  i  a  o 


(6) 


where  the  letters  i,  o  and  a  denote  the  inner  adherend,  the  outer  adherend  and  the 
adhesive  respectively. 

Other  constraints  on  the  displacements  were: 


uj  (O  <  x  <  L,  z  =  0)  =  0 
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u  °\x  =  0,h.+h  <z<h.+h  +h  =0  (8) 

x  '  i  a  i  a  o>  w 

Mesh  refinement  was  carried  out  to  ensure  convergence  of  the  models.  Details  of  the 
mesh  refinement  can  be  found  in  reference  [5]. 

The  FE  value  of  k  was  obtained  by  taking  the  average  of  the  x  displacements  of  inner 
adherend  at  x=0: 


h. 

—  |  u  z(x  =  0  )dz 
ni  0 


(9) 


where  the  integral  was  evaluated  numerically  using  FE  displacements  and  Simpson's 
Rule.  The  average  was  then  inserted  into  equation  1  and  k  determined. 


To  determine  the  range  of  accuracy  of  the  various  analytical  models  as  compared  with 
the  "exact"  FE  results  it  was  first  necessary  to  identify  an  appropriate  measure  of  joint 
configuration  variation.  The  measure  chosen  was  the  non-dimensional  parameter  a 
defined  below: 


«  PpE 


(10) 


Two  length  scales  comprise  a.  The  load  transfer  length,  /?fe,  is  a  finite-element  estimate 
of  the  load  transfer  length.  The  other  important  length  scale  is  the  St  Venant  length 
defined,  for  the  purposes  of  this  study,  as  shown  below: 


St  Venant  Length  =  h 


(11) 


Further  details  on  p  and  the  St  Venant  lengyth  can  be  found  in  reference  5. 
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3.1  FE  Results  for  the  Isotropic  Adherend  Joints 

FE  results  for  a  range  of  isotropic  adherend  DSJs  having  identical  adherends  is 
presented  in  Table  1. 

Table  1.  FE  Values  ofk  and  a  for  the  Isotropic  DSJs. 


h  (mm) 

a 

k  (mnr1) 

1 

0.4818 

0.1303 

2 

0.6618 

0.0887 

3 

0.7853 

0.0698 

4 

0.8822 

0.0585 

5 

0.9625 

0.0507 

6 

1.0289 

0.0449 

7 

1.0862 

0.0405 

8 

1.1358 

0.0368 

9 

1.1780 

0.0339 

10 

1.2127 

0.0314 

7 


DSTO-RR-0139 


3.2  FE  Results  for  the  Orthotropic  Adherend  Joints 

FE  results  for  a  range  of  orthotropic  adherend  DSJs  having  identical  adherends  is 
presented  in  Table  2. 

Table  2.  FE  Values  of  the  k  and  a  for  the  Orthotropic  DSJs. 


h  (mm) 

a 

k  (mm-1) 

0.5 

0.6526 

0.1253 

0.6 

0.7069 

0.1135 

0.7 

0.7564 

0.1042 

0.8 

0.8031 

0.0966 

0.9 

0.8447 

0.0903 

1.0 

0.8828 

0.0851 

2.0 

1.1803 

0.0559 

3.0 

1.3379 

0.0428 

4.0 

1.4708 

0.0351 

5.0 

1.6080 

0.0300 

4.  Analytical  Models 

An  exact  analysis  of  the  stress  state  in  a  DSJ  in  which  both  the  equilibrium  equations 
and  the  compatibility  equations  of  continuum  mechanics  are  satisfied  for  all  three 
layers  and  for  which  continuity  of  tractions  and  displacements  is  enforced  is 
prohibitively  complicated.  Consequently,  simplifying  assumptions  are  made  in 
formulating  an  analytical  approach.  A  trade  off  between  simplicity  and  accuracy  is 
usually  necessary  in  any  approach. 
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The  one-dimensional  analytical  models  of  the  stress  distribution  in  DSJs  considered 
here  fall  broadly  into  two  categories: 

(i)  Those  that  start  from  a  prescribed  form  for  the  displacements  and  in  which 
displacement  potentials  are  introduced.  Strains  are  derived  from  the 
displacements  from  kinematic  relationships  and  automatically  satisfy  the 
compatibility  equations  of  continuum  mechanics.  The  stress-strain  relations  are 
applied  to  determine  the  stresses,  statical  equilibrium  enforced  by  minimising  the 
resulting  potential  energy  and  the  result  is  coupled  ordinary-differential  equations 
for  the  adhesive  stresses.  These  equations  can  be  solved  for  the  adhesive  shear  and 
peel  stresses. 

(ii)  Those  that  start  from  a  prescribed  form  for  the  some  of  the  stress  components  in 
which  stress  potentials  are  introduced.  The  pointwise  equilibrium  equations  of 
continuum  mechanics  are  applied  to  determine  the  remaining  stress  components. 
The  unknown  coefficients  are  determined  by  minimising  the  complementary 
energy  which  ensures  global  compatability.  The  complementary  energy  is 
expressed  in  terms  of  the  stresses  and  the  resulting  functional  is  minimised  to 
yield  the  Euler-Lagrange  equations.  Coupled  differential  equations  are  again  the 
result  and  these  are  solved  for  the  adhesive  shear  and  peel  stresses. 

The  models  described  above  yield  variational  bounds  for  k.  A  further  simplification  of 
the  analysis  results  if  the  condition  of  constrained  elasticity  is  assumed  [6].  Such  an 
assumption  in  this  context  implies  that  for  the  strains  in  the  z  direction: 

s  i  =  s  °  =  s  a  -  0 
z  z  z 

and  furthermore  for  the  z  displacements: 

w  *  =w  °  =w  a  = 0 
z  z  z 

Consequently,  stress  crz  if  it  exists  will  not  contribute  to  the  energy  functional  as  the  sz 
strain  is  zero. 


(12) 


(13) 
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4.1  Models  That  Yield  Variational  Bounds 

4.1.1  Upper  Bounds  -The  Prescribed  Displacement  Models 

4.1.1.1  Hart-Smith 

Hart-Smith's  model  is  a  based  on  simple  plate  theory  and  assumes  that  the 
deformation  state  in  the  adherends  is  one  of  simple  stretching  [7].  The  critical 
assumptions  are: 

1)  That  the  adherends  deform  purely  by  uniaxial  tension  or  compression,  ie,  no 
bending,  so  that  the  only  components  of  displacement  in  the  adherends  are: 

\° W  311(1  uxl(x)  (14) 

2)  The  adhesive  deforms  purely  by  shear  and  the  strain  is  constrained  to  be: 

r M = y  [»x°  M  -  uxl  (*)}  (15) 

a 

The  complementary  energy  of  the  bonded  joint  becomes 


lie  is  to  be  minimised  over  the  range  of  admissible  displacement  functions  subject  to 
the  constraint  given  by  equation  4.1.11.2.  A  Lagrange  multiplier,  X,  can  be  used  and 
the  following  equation  minimised: 


Ua  U  1  ,U  ° ,Y 
11  x  x  '  c 


,y 


a 


y 


a 


(17) 
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This  leads  to  the  set  of  equations: 

A 


d2u  z 


■E.h. 


x  _ 


i 

cl  u 


h  1  1  dx' 
a 


0,--E.h.  . 

*  Z  *  dx2 


X  _ 


=  0, 1+G  h  y  =0. 
a  a'  a 


(18) 


The  first  two  equations  can  be  rearranged  and  subtracted  to  yield: 

h 


( d2u  1  d2u  0 ^ 

JC  X 


dx ' 


dx ' 


a 


f  1  (  1  ^ 

,  E  h  +  E.h.  , 
V  o  o  i  i  J 


(19) 


Substitution  of  equation  19  into  the  last  equation  of  equation  set  18  yields: 


1 d2u  z  d2u 

X  X 


dx' 


dx ' 


a 


-  +  G  h  y  =0. 


/  ^  a  a'  a 

,  E  Y  +  , 

V  o  o  i  ij 


(20) 


Taking  the  second  derivative  of  equation  15  and  substituting  into  equation  20  gives: 


d2y 
■h  a 


a 


a  dx2 


1  1 

+  - 


E  h  E.h. 

V  o  o  i  u 


+  G  h  y  =0. 
a  a  a 


(21) 


This  can  be  rearranged  into  the  form  given  by  Hart-Smith's  [7]  form: 


d*y  G 
a  a 


dx 


2  h 


1  1 
■  +  - 


E  h  E.h. 

V,  o  o  i  i  J 


y  =0. 
a 


(22) 


Thus  minimising  the  complementary  energy  gives  the  same  resulting  differential 
equation  and  as  imposing  statical  equilibrium  (Hart-Smith). 

The  solution  of  which  over  the  semi-infinite  domain  is: 


—  /3x 


(23) 
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where 


G 


a 


a 


E  h 
o  o 


■  + 


i  1 


E.h. 


i  ij 


(24) 


The  constant  yo  is  obtained  by  recognising  that  load  transfer  occurs  purely  by  shear  of 
the  adhesive  and  that  the  load  transferred  equals  cr0h  (in  this  case  of  identical 
adherends): 


JG  y  (x)  =o-Ji , 
J  a'  aK  '  0 

0 


(25) 


Equation  23  becomes: 


,  ^  CT0^  -fit 

a 

k  can  be  obtained  form  equation  \,  recognising  that  at  x=0: 

h 

u  °(o)  =  Oandu  z(o)  =  m  =  — —<y ahj3 , 

X  X  u 

a 

so  that 


E.u  E.h.h  B 

i  i  i  a 


For  identical  adherend  joints  this  reduces  to: 


2 


(26) 


(27) 


(28) 


(29) 


Figure  5  plots  values  of  k  determined  from  equation  28  for  various  joint  configurations, 
confirming  that  the  bound  is  an  upper  bound: 
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Figure  5.  Plot  against  a  of  the  analytical-to-FE  ratio  for  k  -  Hart-Smith. 

4.1. 1.2  Delate,  Erdogan  and  Aydinoglu 

Delale,  Erdogan  and  Aydinoglu' s  [8]  model  of  the  stresses  in  adhesively  bonded  joints 
applies  Reissners  plate  theory  to  the  adherends.  The  adherends  are  allowed  to  undergo 
bending  and  transverse  shear.  Strains  in  the  adhesive  are  assumed  to  be  constant 
across  the  thickness  of  the  bondline. 

The  displacements  ux,  uz  and  the  rotation  f  for  the  two  adherends  are  prescribed  in 
terms  of  traction,  shear  force  and  bending  moment.  Equations  of  statical  equilibrium 
relate  the  traction,  shear  force  and  bending  moment  to  the  adhesive  shear  and  peel 
stress. 

Delale  et  al  derived  the  following  coupled  differential  equations: 

3 

dr  dr 

- -£-+<x* — —+a~cr  =  0  (30) 

1  dx  2  a  > 


4  2  3 

d  o  d  o  d  r 

(31) 
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The  constants  are: 


a,  =~ 


G„ 


h 


a  \ 


C1+C2  +  (ho+h°)h°D  ^ 
4 


a2  =— D0, 

K  2 

[2(1 

D0 

* _ _  VA 

2(l-v.-2vS 

R  —  ' 

2Aa(l-va-2v„2; 

'jL+i^ 

vA- 


C1-C2  +  ^^0+^)i)o 

4  ; 


where. 


r  -1z2Zl  r  -l~viv‘ 

°  K  E,  ’  1  hft  ’ 

_  12(l-iyv,) 

0  El 

The  terms  in  those  constants  that  relate  to  bending  of  the  inner  adherend  (D;=0)  can  be 
set  to  zero  for  the  DSJ.  The  procedure  for  solving  these  coupled  ordinary  differential 
equations  with  constant  coefficients  is  straightforward:  reducing  the  two  equations  to 
one  higher  order  equation  in  one  of  the  stresses  and  substituting  an  exponential  of 
form  ex*  and  solving  the  resulting  characteristic  equation  for  X.  However,  the  roots  of 
the  characteristic  equations  are  quite  complicated  and  messy.  By  numerically 
evaluating  the  constants  and  using  the  program  Mathematica  [9]  exact  solutions  can  be 
found  though  not  in  symbolic  form. 
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Values  of  k  for  the  various  joints  are  shown  in  Figure  6,  also  indicating  the  expected 
upper  bound: 


Figure  6.  Plot  against  a  of  the  analytical-to-FE  ratio  for  k  -  Delate  et  al. 

4.1.2  Lower  Bounds  -  The  Prescribed  Stress  Models 
4.1. 2.1  Hashin 

The  analysis  developed  by  Hashin  [3]  for  cracked  composite  laminates  is  extended  to 
bonded  joints  in  this  section.  Admissible  stress  states  for  the  adherends  and  adhesive, 
satisfying  the  equilibrium  equations  of  continuum  mechanics,  are  developed  in  terms 
of  a  single  stress  potential,  (f>.  The  complementary  energy  of  the  joint  is  minimised  to 
yield  an  ordinary  differential  equation  in  <f> .  Solution  of  the  differential  equation  allows 
the  determination  of  an  expression  for  k. 


The  assumed  stress  state  is: 
Inner  Adherend: 


= -<?Ax\ 

=  M'kk  (32) 

=  cr^\x)\  (hi  (hi  +  ho)-z2} 
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Outer  Adherend: 


Adhesive: 


<JZX°  =  croT~ 0'(xXhi  +  K  +K~  z\ 

no 

+K+K  -zf- 

=  o, 

vj  = 


The  complementary  energy  is  given  by: 


hi+h0+ha 


Uc-  ^dz^W(x,z)dx, 


2W  =  ^(^  +  ^  +  ^--2fcrxxcrzz  . 

V  Et  Et  G  t[  J 


Contributions  from  the  adherends  and  the  adhesive  are  summed.  This  yields: 

2  Uc  =  <J2  h2  jdx{d(,^2  +  A  (&'Y  +  AH>"  +  A  iff  }> 

0 

where, 

AizkiiL  ,l~Wl 

4  Et\  e:\  ’ 

A  =  ik+ A_+A_ 

1  Ga  3  G,l  3  Gu° 

^ A'++,(l+,/'')'(i  v)^/(l+v;) 

.  (  i  V  \  ^  .•  N  1.3(1  (. .  0  V  ^  1  l  (1  , ,  2  \ 


4 E!  1 15  3  ‘  0  )  20 E  4  Ea 


Minimising  Uc  leads  to  the  following  differential  equation: 


At""  +  (A~AW  +  A  =  o- 
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The  constants  f3i  and  [k  are  given  by: 


From  this  solution  for  <|)(x)  the  stress  state  for  the  joint  is  determined,  k  can  be  found  by 
integrating  the  complementary  energy  (equation  4.3.?): 

'  2  A&fe+A) 

k  is  then  found  from  equation  2. 
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Values  of  k  for  the  various  joints  plotted  in  Figure  7  below. 


Figure  7.  Plot  against  a  of  the  analytical-to-FE  ratio  for  k  -  Hashin. 


4.1.2.2  Allman 


Allman's  analysis  [10]  of  the  double-lap  joint  utilises  two  stress  potentials, 
<f>i(x)  and  to  describe  the  stress  state  of  the  joint,  viz: 


1 

2 


2*  ex 2 


f 

2 


r 

<7 

o 

\ 


\ 

+  <7  . 

I 


(43) 


(44) 


where  cr0  and  cn  are  the  stresses  normal  to  the  bondline  at  the  interface  with  the  outer 
and  inner  adherend  respectively. 

Allman  models  the  effects  of  bending,  stretching  and  shearing  of  the  adherends  and  of 
shearing  and  tearing  of  the  adhesive  layer  using  these  two  stress  potentials.  The 
advantage  of  two  stress  potentials  over  one  is  that  requirement  of  a  stress-free 
condition  at  the  end  of  the  adhesive  layer  can  be  accommodated: 


<C(o)=o. 


(45) 
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The  form  of  the  stresses  in  terms  of  stress  potentials  that  Allman  derives  [10],  with 
some  modification  to  allow  for  the  infinite  nature  of  joint,  for  the  double-lap  joint  is: 


for  the  outer  adherend 


of  \  1  2 


K  (K+hi) 


(3(1 +£,}]- 1 V,  M+  6r}f2^ , 


r-”w=iM((1+3£')+2,?'3(1+',>72^+^:(l"'?2)^’ 

ayy°  M  =  °+  h^j  (0  +  2s1)-{l  +  3sl  yj-rj1  +(l  + 

If,  3  1  3VV2 

2{  2  2  J  dx2 

where  ij  =  ho/2  and  £i=h/h0 


for  the  inner  adherend 


where  rj  =  h/2 


1  4 


h,K+h/' 
dfa  7J+ 1 
**  dx  hi+hg’ 


'  dx2  dx2  *i+A0 1°  4l  7  V  >)’ 


for  the  adhesive 


«  2 
^  dx  h,  +  h  ’ 

a ‘=£A+„J±-££.  (48> 

w  dx2  7  \  +ho  dx1  ’ 

where  77  =  h/2. 

The  strain  energy  of  the  joint  is  expressed  in  terms  of  the  stress  potentials  and  then 
integrated  over  z  and  x.  The  resulting  functional  is: 
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an  ~ 

0C13  = 


i  (l+e^Xi-VoK 

ftr+O  5  E0 

6(1- 

s*A 


which  is  rniriirnised  to  yield  Euler-Lagrange  equations  for  the  two  potentials.  These 
coupled  diffemtial  equations  are  of  the  form: 


2 («,  +  (a10  - a4)D2  +  a7D 4)^  +  ( a3  +  {an  +  al2  -  a6)D 2  +  a9D*)<j>2  =  0  (50) 

and, 

(«3  +  («n  +  an  - a6)D 2  +  a9£>4)$  +  2(a2  +  (an  - ccs)D2  +  a8D4)^2  =0  (51) 


These  equations  can  be  combined  to  give  the  following  eigth-order  differential 
equation  for  <j)2: 


d\ 


y  d\  d%  d\ 

1  2  2  dx2  3  dx4  4  dx6  5  dxS 


where  /?/  to  J35  are  related  to  a\  to  an  by: 
=  -4«,a2  +  a2 


(52) 


P2  =4(«,(«5  -a13)+«2(«4 -«io)+^3(«n  +«12  -««)) 

^3  —  1 1  "i"  ^12  —  ^6  )  —  4((o,0  —  X®13  —  ^5  )  ^2*^7 

y?4  =  -4a7  (al3  -  a5 ) +  4a8  (a10  -  a4 ) +  2a9  (a, ,  +  a12  -  a6 ) 

=  -4a7ag  +  a92 

Equation  52  involves  only  derivatives  of  even  powers  and  so,  in  effect,  reduces  to  a 
fourth  order  differential  equation.  Consequently,  solutions  for  <|>2[x]  that  decay  to  zero 
as  x  approached  infinity  are  of  the  form: 
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<t>,  [x]  =  AE~X^+  BE~X &  +  CE~X'^  +  DE~X^  (53) 

where  the  constants  X,  to  X4  are  real  roots  of  the  quartic: 

+  /32u  +  /?3«2  +  /?4u3  +  y 3su4  =  0  (54) 

■The  roots  X,toX4  all  had  positive  real  parts  for  the  joint  parameters  investigated. 
Typically,  Xt  to  X4  are  real  and  X4  and  X4  are  complex  conjugates.  The  stress 
potential,  <j>i[x],  has  the  same  form  as  <|>2[x]: 

A  [x]  =  EE~X^  +  FE~X &  +  GE~X ^  +  HE'X &  (55) 

The  constants,  E,  F,  G  and  H  are  related  to  A,  B,  C  and  D  by: 


J7  .  .  3  +  ( au  +  a\2  )^1  + 

2(«i  +  (^io —  )2j  +  ct-jX^ ) 

(56) 

F  Ba^Jr  +  an  ~  + 

2(al+(axo-a4)X2+a1X2) 

(57) 

(J_Cai+  11  ®12  ~  ^6  )^3 

2(al+(al0-aJX3+a7Xi) 

(58) 

JLJ  _  £)a3  +  (®11  +  ®12  ~  a6  )24  +  <^9^4 

2(«i  +(aio-«4K+M4) 

(59) 

The  boundary  conditions  applied  to  the  two  stress  potentials  were  that: 


r/(0)-0^A[0]  =  0 

(60) 

Va(x) 0  as  x  -* oo  =>  ^^-[0]  =  0 

dx 

(61) 

K'  (x)dx  =  -F  =*[0]  =  -F 

0  ^ 

(62) 

M„(x)-> Oasx  ->  oo  =>  d<^2  [0]  =  0 

dx 

(63) 
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where  M0  and  V0  are  the  bending  moment  and  shear  force  resultants  in  the  outer 
adherend  respectively. 

Once  the  two  stress  potentials  are  determined  k  can  be  found  by  integrating  Uc  directly 
via  equation  2. 

The  results  are  shown  in  Figure  8  for  the  two  types  of  joints. 


Figure  8.  Plot  against  a  of  the  analytical-to-FE  ratio  for  k  -  Allman. 


4.2  Constrained  Elasticity  Models 

The  constrained  elasticity  models  impose  the  additional  constraint  (transverse- 
inextensibilty)  that: 


£z=0  .  (64) 

Constrained  elasticity  models  have  been  used  successfully  by  Savoia  and  Tullini  [11]  to 
derive  a  beam  theory  for  strongly  orthotropic  materials.  This  approach  can  be  used 
when  the  height  of  the  beam  is  sufficiently  small  with  repsect  to  the  length  or  when  the 
Young's  modulus.  The  stress  tensor  then  divides  into  an  active  and  a  reactive  part.  The 
reactive  part,  <rzz ,  does  no  work  for  any  admissible  deformation. 
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For  bonded  joints,  which  are  typically  long  and  slender  and  often  have  orthotropic 
adherends,  transverse  inextensibility  applies  to  both  adherends: 

u',0z,z  =  0=>  u‘'°z(x,z)  =  v‘‘0(x).  (65) 

Symmetry  imposes  the  following  boundary  condition: 

u'z  (x,0)  =  v1  (x)  =  0 .  (66) 

Within  the  constrained  elasticity  format  it  is  possible  to  prescribe  the  format  of  either 
the  displacements  or  die  stresses.  The  guarantee  of  true  upper  or  lower  bounds  is  lost 
however.  These  two  approaches  are  described  below. 

4.2.1  Prescribed  Displacement  Model 

The  form  of  the  displacements  in  the  adherends  can  be  specified  using  the  polynomial 
form  adopted  by  Savoia  and  Tullini  [11]: 

«'.(*.<r)=2XfeK,( 4-  (67) 

o 

where  £  =  2z/h,. 

The  polynomials  Pm  are  defined  by  a  recurrence  relation: 

■P.-l.  (68) 

Pn  also  has  the  useful  property  that: 

i  i 

[P.d{=f&,df  =  0forn>2.  (69) 

-1  -1 

This  leads  to  the  interpretation  of  ®o(%)  as  the  thickness-averaged  value  of  u{x. 

A  solution  to  order  P2  can  be  obtained  as  follows.  For  the  inner  adherend  the 
displacements  are: 

u‘x  (x,  %)  =  ui  (x) + Px  (x) + P^yp‘2  (x) 

=  u‘  (x)  +  £0'i  (*)  +  <J>'2  (x) ,  (70) 

6 

u'z(x,%)  =  0. 
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The  shear  stress  in  the  adherend  is  then: 


<*' *(*,4)=  (*)+  i>'(f)l>'2(xj 


•The  boundary  condition  on  the  inner  adherend  shear  stress  is: 


which  implies. 


So  that. 


ar'xt  (x,  £  =  ±l)  =  ±r(x). 


®!2(x)  =  — LT(x)and®i(x)  =  0. 
G, 


u‘x (x,  £)  —  u‘  (x)  -  ^  * t(x)  , 

(j;  O 


Similarly,  for  the  outer  adherend: 

U°x  (x,  £)  =  u°  (x)  +  fe)®0,  (x)  +  P2 (<^)002  (x) 

=  u°  (x)  +  %$>\  (x) + ^ ■— ^  0°2  (x)  ' 

6 

The  boundary  condition  on  the  outer  adherend  shear  stress  is: 


So  that. 


cr'xf  (x,  £  =  ±l)  =  +r(x) . 


O02  (x)  =  — —  r(x)  and  (x)  =  0, 

G0 

u°x{x,^)=u°(x)+^-  r  (x) 


As  for  the  Hart-Smith  analysis,  the  complementary  energy  of  the  bonded  joint  becomes 
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(78) 


■  lie  is  to  be  minimised  over  die  range  of  admissible  displacement  functions  subject  to 
the  constraint  given  by  equation 


y(x)=  ~(ux°(x,%  =  =  l)) 


a 


t  (  h.  h„ 


W 


=  —  u  -u - — -  + 

3  G, 


'oJJ 


(79) 


A  Lagrange  multiplier,  X,  can  be  used  and  the  following  equation  minimised: 


U,  (O u° ,  ra  )=  Uc  (ux‘  ,ux°  ,ya )+  A 


Ya  — 


h  h  h 

— - — | H 2_ 

'\Gm  G,  Ga 


■(u°-u‘) 


(80) 


This  leads  to  the  set  of  equations: 


X  7  dlu  '  X  d2uj 

T  -  M  -Ti-  =  0  -  -  EA  -TT  =  0-  ^GahJa  =  0. 


h  1  ‘  dx2  ’  h 


dx2 


(81) 
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These  equations  can  be  solved  to  yield  solutions  for  the  displacements,  k  can  be 
determined  from  the  same  approach  as  used  in  the  Hart-Smith  analysis. 

Figure  9  plots  values  of  k  determined  from  equation  2  for  various  joint  configurations: 


4.2.2  Prescribed  Stress  Models 

Prescribed  stress  models  can  also  be  developed  within  the  constrained  elasticity 
framework.  Two  models  are  developed  -  one  that  models  adherends  tractions  only  and 
one  that  models  tractions  and  moments. 

A  slighlty  different  nomenclature  is  used  in  the  following  two  models:  the  inner 
adherend  thickness  is  taken  as  hi  rather  than  2/i,  as  previously. 

4.2.2. 1  Tractions  Model 

The  adherend  stresses  are  defined  in  terms  of  inner  adherend  traction,  Ni,  and  the 
outer  adherend  traction,  N0: 


<t‘x(x)  = 
(7°x(x)  = 


Nfc) 

hi  ’ 

K(*) 

K 


(82) 
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The  other  active  stresses  are  determined  from  the  equilibrium  equations: 


aixz(x,zi)  =  -2zi^^-, 
cr0X2{x, za)  =  N0  (x0-^ 


(83) 


where  z,  and  z0  are  coordinate  systems  local  to  the  inner  and  outer  adherends 
respectively. 


The  boundary  conditions  are: 

crlxz(x,  ^  /  2)  =  - N ]  (x), 
a‘xz(x,0)  =  0, 
(7°X2(x-h0/2)  =  N^(x), 

a°xz(x,+h0/  2)=0, 


(84) 


The  complementary  energy,  equation  3,  is  constructed  and  minimised  and  the  resulting 
Euler-Lagrange  equations  solved  to  give: 

ya(x)=roe~lk  /  (85) 
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where,  for  an  orthotropic  material: 


^(l-tW) 

(i-V/XT 

( K  +  K 

+  K  ^ 

I 

E  °  h 

V  no 

E/h,  , 

3  Gt°) 

(86) 


k  is  given  by  equation  29  for  an  identical  adherend  joint.  Plots  of  k  are  shown  in  Figure 
10  below. 


4.2. 2.2  Tractions  And  Moments  Model 

The  adherend  stresses  can  be  defined  in  terms  of  tractions  and  moments: 
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The  other  active  stresses  are  determined  from  the  equilibrium  equations: 

N'(x) 


(Tixz{x,Zi)=-2zi  ■ 


h 


cr°xz(x,z0)  =  Na  {x\\-?A+M0  (*)V 


2  h 


O  J 


2  L 


(88) 


The  boundary  conditions,  equations  84,  are  satisfied  by  these  equations. 

The  complementary  energy,  equation  3,  is  minimised  and  the  resulting  Euler-Lagrange 
equations  solved.  The  resulting  expressions  were  very  complicated.  By  numerically 
evaluating  the  constants  and  using  the  program  Mathematica  exact  solutions  were 
found  though  not  in  symbolic  form. 


k  can  be  determined  from  integrating  the  complementary  energy.  Plots  of  k  are  shown 
in  Figure  11  below. 


Figure  1 1.  Plot  against  a  of  the  analytical-to-FE  ratio  for  k  -  S&V/Rose. 
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5.  Discussion  and  Conclusions 


Expressions  that  yield  variational  bounds  for  the  equivalent  spring  constant  of  the 
double-strap  joint  have  been  obtained.  These  expressions  should  be  useful  in  obtaining 
more  exact  estimates  of  the  stress  intensity  factor  associated  with  patched  cracks. 

True  upper  bounds  were  obtained  from  the  Hart-Smith  and  Delale  et  al  displacement 
based  models.  The  Hart-Smith  model  is  used  in  RAAF  Engineeering  Standard  C5033 
and  gave  the  least  conservative  results.  Delale's  model  gives  improved  results  but  is  a 
mathematically  complicated  formulation. 

True  lower  bounds  for  the  equivalent  spring  constant,  which  are  useful  for  obtaining 
truly  conservative  estimates  for  patched  cracks,  were  obtained  from  the  Allman  and 
Hashin  stress-based  models.  The  Hashin  model  was  more  accurate  for  the  joints 
considered,  especially  for  orthotropic  joints. 

The  constrained  elasticity  models  worked  quite  well  although  the  guarantee  of  true 
upper  or  lower  bounds  was  lost.  However,  the  tractions  and  moments  model  gave 
conservative  results  over  the  range  of  joints  considered. 

The  Hashin  model  was  the  probably  the  best  model  considered  -  it  gave  truly 
conservative  estimates  of  k,  was  reasonably  accurate  over  a  wide  range  of  joint 
parameters  and  an  analytical  expression  was  obtained  for  k. 
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